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(57) Abstract 

A method and apparatus are provided for compressing data 
received from a digitizer tablet based on the characteristics of each 
individual stroke contour (150, 160). Tlie digitizer tablet samples 
the position of the writing pen, continuously transmitting data to a 
computer in the form of x and y coordinates plus an indicator of 
whether or not the pen is touching the surface of the tablet (100, 
1 10). After preprocessing, a determination is made of which points 
along the stroke contour are essential to the reconstruction of a 
quahty facsimile of the original writing from the sampled points 
(140). This determination is based upon local cur\ature, local ex- 
trema, and the endpoinis of the stroke. In an optional second stage 
of compression, a standard mathematical compression technique 
is applied (160). 
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BACKGROUND OF THE INVENTION 
The present invention relates generally to data 

capture applications, and more specif icallv tn i-y,^ 

r,-F ^ = f«=<-Aiiv,axiy, to the compression 

of data generated on a digitizer tablet. 

10 When something is written on a digitizer tablet the 

output zs a parametric representation of the writing 
trajectory; that is, the writing is represented as a series of 
X, y coordinate values as a function of time, a typical 
dxg.ti,er tablet samples the position of the writing pen loo 

.IT/" ''''' '"^^^^^^ °^ -"n^- The 

.ablet continuously transmits the resulting data to the 

la'mo ifL'"'"'"' "'^^ ^^^'^ tablet 
samples the pen position and transmits an x-coordinate value a 

Whether the pen is up or down. The data generated in this 
manner is commonly referred to as electronic ink. 

There are applications where it is of interest to 

oes^raMe t= save an electronicaUy generated signature Lr 
!!! -"^niticn. Another application is the 

P r^anent storage of aata whioh has been generated on a 
digitizer tablet. 

s^orage J^.re^ or.e^X rgYT/ thf '' '''' ""^ 

J ^ -ue erricient to compute. The 
c=.,p.tarional efficiency is important because electronic ink 
Zl °" ''^'"''"^ "^"^-^y -^-t CPU pc«r 

: °' '^"'^"^ Significantly Slowed if 

.Pie. aeco.pression techni^e had to be applied to each ^en 

-cTbegLTi: rpftiht"::: t 

snr-^r-o . - . ^oucnes the digitizer tablet's 

surface ana ends when the pen is * 

pen is lifted from the surface ) 
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There are a number of methods for data compression. 
The earliest methods were based on the idea of minimum 
redundancy coding. The basic idea was very simple; symbols 
which occur frequently are coded with fewer bits than those 
which occur more rarely. The earliest known method to apply 
this technique was the Shannon-Fano coding. A subsequent 
development was Huffman coding, which is similar to Shannon- 
Fano coding in many ways. Huffman's contribution was in 
devising a very clever and elegant way to construct the coding 
and decoding trees. In general, simple Huffman coding does not 
achieve the same degree of compression as more modern 
techniques, but processing and memory requirements are 
relatively modest. 

A further enhancement of Huffman coding is called 
adaptive Huffman coding. A basic problem with standard Huffman 
coding is that the decoding program must have access to the 
Huffman coding statistics. Hence, trying to increase 
compression with more comprehensive statistics is eventually 
self-defeating since more and more coding statistics must also 
be saved, which works against the compression. Adaptive 
Huffman coding gets around this problem by adjusting the 
Huffman decoding tree on the fly, effectively gaining the 
advantages of higher order modeling with -3 additional storage 
for statistics. 

Huffman coding requires that an integral number of 
birs be used to code symbols. Arithmetic coding improves 
compression by using a fractional number of bits per code. 
Although arithmetic coding will always provide as good or 
better compression than Huffman coding, its disadvantage is a 
significantly larger computational burden. 

Dictionary coding is another class of data 
compression. In general, these techniques create a dictionary 
of symbol sequences or combinations which are then accessed by 
indexing into the dictionary. Th'ese techniques have become the 
sranc -d for compression on small computers because they 
combi : good compression with modest memory requirements, 
Examp.es of codes using these techniques which are applicable 
to MS-DOS include PKZIP, ARC, ARJ , and Lharc. 
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T^^^® also compression techniques that apply 

primarily to specialized domains. An example is compression of 
two dimensional images using the discrete cosine transform and 
its many derivatives. 

A method and apparatus for compressing electronic ink 
which achieves high compression while requiring minimal memory 
and CPU power is needed. 
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SUMMARY OF THE INVENTION 
According to the invention, for a computer with a 
digitizing tablet, a method and apparatus are provided for 
compressing data based on the characteristics of each 
individual stroke contour. 

The digitizer tablet samples the position of the 
writing pen a pre-determined number of times per second and at 
fixed intervals. This data is continuously transmitted to a 
computer in the form of x and y coordinates plus an indicator 
of whether or not the pen is touching the surface of the 
tablet. 

In the preferred *jmbodiment the incoming data from 
the digitizer tablet is preprocessed. This preprocessing, 
which is application and/or system specific, is not part of the 
compression. Examples of preprocessing include smoothing for 
noise reduction, andpoint filtering, discarding duplicate data 
points, median filtering. Manning filtering, spatial sampling, 
and size noraalization. 

In the first stage of compression a determination is 
made of which points along the stroke contour are essential to 
the reconsrrucrion of a quality facsimile of the original 
writing from the sampled points. The method proceeds by 
calculating descriptors and finding essential features of the 
stroke contour. In the preferred embodiment, the descriptors 
are local curvature, local extrema , and the endpoints of the 
stroke. Although the actual algorithm is somewhat detailed, in 
general the idea is that the extrema are saved, the points ' 
about Which the local curvature is above a threshold are saved, 
and the endpoinrs of the stroke are saved. 
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^" second stage of compression any number of 
standard compression techniques can be applied including 
Huffman coding, adaptive Huffman coding, arithmetic coding, and 
a variety of dictionary compression techniques. In the 
preferred embodiment differences of the form Ax = - x^.^^ and 

= Xi - Vi-i are calculated and the stage 2 compression is 
applied to the differences. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a flow chart of the general compression 
method according to the invention. 

Fig. 2 is a flow chart of the extrema determination 

step. 

Fig. 3 is a flow chart of the stroke endpoint 
determination step. 

Fig. 4 illustrates the local curvature calculations. 
Fig. 5 is a flow chart of the duplicate data points 
preprocessing step. 

Fig. 6 illustrates the endpoint filtering 
preprocessing step. 

Fig. 7 is a flow chart of the endpoint filtering 
preprocessing step. 

Fig. 8 is a flow chart of the median filtering 
preprocessing step. 

Fig. 9 is a flow chart of the Hanning filtering 
preproceissing step. 

Figs. lOA-C illustrate the reconstruction of an 
electronically written signature without applying this 
invention's compression techniques. 

Figs. iiA-c illustrate the reconstruction of an 
electronically written signature after applying this 
invention's compression techniques. 

Figs. 12A-C illustrate the reconstruction of a 
sequence of lower case letters after applying this invention's 
compression technique. 

Figs. 13A-C il-ustrate the reconstruction of a simple 
line after applying this i.nvention's compression technique. 
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Figs. 14A-E illustrate the trade off between 
compression ratio and reconstruction fidelity. 

Fig. 15 illustrate the apparatus to accomplish the 
disclosed method. 

DESCRIPTION OF THE SPECIFIC EMBODIMENT (S ) 
The flow chart of Fig. l illustrates the compression 
method described by this invention. Data is input by writing 
on a digitizer tablet (step 100). The tablet represents the 
writing as a series of x and y coordinate values as a function 
of time (step iio) . This stream of x,y position coordinates 
along with a pen up/pen down indicator for each coordinate pair 
is continuously sent to the computer (step 120) . 

In the preferred embodiment, prior to data 
compression, the data is preprocsssed (step 130) . The exact 
form of preprocessing will vary according to the application 
and system. Preprocessing is generally used to remove noise. 
The details are unimportant as long as a good representation of 
the written input data is maintained. Another optional step 
which is dependent upon the application and system is the 
removal of duplicate or close points (step 140). This step is 
generally necessary when the digitizer tablet's user writes 
slowly. Slow writing causes oversampling resulting in duplicate 
or close points which are redundant and do not carry 
significant information. 

The first stage of compression is based on the 
characteristics of each individual stroke contour (step 150). 
In the second stage of compression any of the standard 
mathematical compression techniques can be applied (step 160). 
Although this stage of compression is not required by the 
invention, it is employed in the preferred embodiment as a 
potential means of further compressing the data. 

After the completion of the second stage, a 
comparison is made between the first stage and the second stage 
compression (step 170). if the compression was not improved 
after applying stage two, then the stage one compression data 
IS saved (step 180) . if the data was further compressed by 
stage two, then this compression data is saved (step 190). 



^'^^''/OaSSa PCr/US93/06883 

6 

The computer program modules used to compress the 
electronic data are given in the appendix. Significant program 
modules and their specific functions are as follows: 

Stage one compression: 

caic.sav.c Function to determine which points in the 

signature need to be saved, based on various 
criteria. 

curvature. c Function to calculate the local curvature at a 

point along the contour. 
Stage two compression: 

cmpsig2.c Function to compress an array of stage 1 data 

with characteristics of the signature 
verification template data. 

pack_hal f bytes, c 

Function to pack the value to be compressed into 
half bytes. 

set_4_bits Function to pack values into the upper or 

lower 4 bits of a byte. 
get_4_bits Function to unpack values from the upper or 

lower 4 bits of a byte. 
uncmpsig2.c Function to uncompress an array of data 

compressed by the function cmpsig2.c. 
unpack_half bytes . c 

Function to unpack and reconstruct the value 
from half bytes. 

Figs. 2-4 are detailed illustrations of the stage one 
compression methods. Fig. 2 shows a flow chart of the first 
step Of compression in which extrema are determined and saved. 
The axg.t.zer data is in the form of a sequence of coordinate 

aiues, Xi and y., as a function of time (step 210). The 
subscript i varies from i=o to i=nsamps-2 where nsamps is the 
total number of data points in a stroke. The subroutine is 
initialized by setting i equal to l (step 220). 

Minima and maxima are determined by comparing- the 
value Of each data point to the values of the data points 
directly adjacent to it (step 230). To be a maxima, the value 
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"of either the x coordinate or the y coordinate must be larger 
then the values of the corresponding coordinate of the data 
points on both sides of the data point under review (step 230) . 
To be a minima, the value of either the x coordinate or the y 
coordinate must be smaller then the values of the corresponding 
coordinate of the data points on both sides of the data point 
under review (step 230) . 

If any of the four relationships defined by step 230 
are met (step 240), then that data point will be saved (step 
250) , otherwise the data point is not saved. After determining 
whether or not to save a data point, the value of i is 
increased by one (step 260) and a determination is made as to 
whether or not all data points have been reviewed (step 270). 
If all data points have been reviewed then this subroutine ends 
(srep 280), otherwise the subroutine loops back to the 
comparison step (step 230), continuing until all data points 
have been reviewed. 

Fig. 3 shows a flow chart of the second step of 
compression in which stroke endpoints are determined and saved. 
The digitizer data is in the form of a sequence of coordinate 
values, xj_ and y^, as a function of time (step 310). The 
subscript i varies from i=0 to i=nsamps-l where nsamps is the 
total number of data points in a stroke. The subroutine is 
initialized by setting i equal to 0 (step 315). 

This subroutine first determines whether or not the x 
and y coordinates of the data point indicate that it is the 
first data point of the stroke (step 320) . If it is the first 
data point, then it is saved (step 330) . Next the subroutine 
determines whether or not the x and y coordinates of the data 
point indicate that it is the last data point of the stroke 
(step 340). If it is the last data point of the stroke, then 
it is saved (step 350) . After determining whether or not to 
save the data point, a comparison is made to determine if all 
data points have been tested (step 360) . If they have been 
then this subroutine ends (step 370), otherwise the value of i 
is increased by one (step 380) and endpoint determination 
continues . 
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^^9- ^ illustrates the calculations used to determine 
the curvature at each data point between the stroke endpoints. 
The stroke contour is represented by line 410 while the 
individual digitizer data points are represented by points 420. 
Point 430 is data point i, the data point of interest in this 
illustration. The coordinates of 430 are and y-. 

In Fig. 4 line 440 is defined as the line segment 
connecting point 430 with point 425. The coordinates of point 
425 are x^.^ and y^.^, where d is equivalent to 3. Line 450 is 
defined as the line segment connecting point 430 to point 435 
where point 435 has coordinates Xi^^ and y^^^. curvature angle 
460 IS the angle defined by lines 440 and 450. 

During stage one compression, the curvature angle for 
each data point is calculated. If the calculated curvature 
angle is within a certain range of angles, the data point is 
saved. If the angle does not fall within the pre-selected 
range of angles, the data point is discarded. The choice of 
the variable "d" as well as the range of threshold angles is 
determined by the application's requirements, m the preferred 
embodiment, the value of d is 2 ; data points which have a 
curvature angle between 5 degrees and 25 degrees are retained 
for applications requiring high fidelity, low compression 
(e.g., signature verification); and data points which have a 
curvature angle between 25 and 50 degrees are retained for 
applications requiring low fidelity, high compression (e.g., 
simple data storage) . 

Figs. 5-8 are detailed illustrations of optional data 
preprocessing methods. Fig. 5 is a flow chart of the method by 
Which duplicate data points are discarded. The digitizer data 
IS m the form of a sequence of coordinate values, x^ and y- 
as a function of time (step 510) . The subscript i varies fJom 
1-0 to i-nsamps-l where nsamps is the total number of data 
pomrs in a stroke. The subroutine is initialized by setting < 
equal to 1, x^^^ equal to Xq and y^^^ equal to y^ (step 520). 

The current data point, represented by x^ and y- is 
compared to the previous data point, represented by x^^^^and 
Vsav (Step 530). If the values for the x and y coordinates are 
idenrical, then the current point is discarded (step 540) and 
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'^the next point is compared to x anri 
continues either aU o. l^MT. ! ! Zr 

reviewed (step seo, or until L no etai . an. 

ysav. respectively. «hen the t« data points ar. f^" 

equivalent, x 



" X33^ is set equal to x. and y is s.t / T . 

(step 545). The value of i is tho^ • ^° 

;,r,^ increased by one fsten «;sn^ 
and the new value af i tstep 550) 

points (Step : . T^e su3" " ""-^^ 

points have\een'.evL.raT rL^tl" T" "^"^ 

10 (step 570). ^® subroutine ends 

Often the beginning and ending of a stroke hav. 
spur.ous artifacts which can be removed with little .... 
the writing's fidelity. These artifa../ t °" 
writer's unsureness Lul. ""^^^"^^ '^^^ to the 

^3 contact, or duerpe:i\" ^Tl^T 

preprocessing .ethod of endpo.!; fill' ^''"^""^^^ 

remove these artifacts. ^^^ter.ng which can be used to 

Stroke contour eio 
points 620. Circle 630 °' " '^"^ 

0 point 635. and clicie lo r " ""^ """^'^ ^"^ 
point 645, are used to aei' " "^"^ 

either discarded slT VlV'^T' 

Physical radius of this circl! ^"^^^^^^ e„bodi„ent, the 

nunber of data oointf I " """»"ers. The 

I data points contained within these circle. i= 
' dependent upon the resolution of the digitizer taWet ™ 
resolution of a typical ■ ^^itizer tablet. The 

Therefore the threshold t " »iUi.eter. 

nuitipiyin, the de L 'ph":::.:",::,"" " """-^■^ 

the tahlet.s resolution XiL^^ute^e^r: 

Preproces/^l- Jethod^. ^e r^ti^r 1^^"^"""^' 
seouence of coordinate values " aL v °' ' 

<«ep ,10,. The subscript < ;aries 1 " ■ " ' °' ""^ 

"here nsarnps is the tot.i " I "° '° i-"5=»Ps-l 
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-indicating that this point is within circle 630, the point 

step 73 ^^'^^ ^"'^ calculation for d 

(Step 730) IS repeated. This process continues untn h • 
5 less then the threshold value at whi.h T !! 

Yj is saved as the second da^I °' ' 

J , second data point (step 770) . This 

(Step 715) and is set equal to nsainps-2 (step 720^ 

su.==:;p:\ i; " ""^ '"-^ • 

" t=«l „u,*er =f data ^""""^-^ ""er. nsaaps is the 

initialized by s«tinr ° ""^ -^-"""e is 

^ setting i equal to l (step 820) 

=ocrdinate"is' «irute7r'T ' 

successive data point °" °^ '""e 

^= tnen in reas d , I ' ^^^^ ^ 

until the V lue o^ i ' '"^ 

-inus two (s ep slo H Tl """"^^ " 

(Step 870, . '^"^ subroutine is complete 

:-"-^^^":^rren;^rirtre:ry^ - 

subsc ;pt\ var i; r ' °' '"^'^ • 

total nuL ; r , rii::: " ™''-' --^^^ "^-^^ ^'■^ 

3. initialled b. J^^^TZZ r/Hsto Z'''''"''- 

-d y coordinate is calculated based on the initial values o. 
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"three successive data points (step 930 ) Th. , , 
then increased by one (step 940) Ld t^a . °' ' 

until the value of i i. » ! subroutine continues 

-nus two (ste ll ZTV''^: """^^ " 

5 (step 960). ^^""^ ''''^ subroutine is complete 

°f storage. The points UllT 

reconstruction fidelity. Fig lie is "^"^ " "^"^ 

signature based upon the cir!i.. -instructed 
required only =1. byt^s of "t ™" 
original. The stage !„ " " °^ 

- ^nes, „hich":: :::t::~: tTr'" - 

«age two compression technique. " ''^ ""P^^^"' ^ 

Fig. 12A is the original nf = — 
letters. As in Fig. ub Fio I «1"^nce of lower case 

1210 received f-o. the . ' """^^ ^=1"= 

- determined to Litlcal t"":"' ^"^ -^"^ 

-c is the reconstruct": seTu::::::::::"" "''""^ 

original required 2072 bytes the rL ! 
-.^ one compression only r:,:^™"::^:,""^ 
the original) , and the reconstruction k °^ 
» both stage one and stage two =clr' ■ " 

201 bytes ,15 percent) «« applied required 

auto.aticany adTul'": te"" T" 

> little structure Fi„ =""9ht line which contains 

"cture. Fig. 13b shows both aci-n»i ■ • 
received from the digitizing t.h, T ^ 

to be critical to th! reco„\ ""'"^ determined 

to tne. reconstruction fideiii-rr r. • 
reconstruction based upo-: circled ".^^^^^"^y* ^^5' is the 

PC circlea points 1320 in which only 



Fig. 
The 
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'■ stage one compression has been applied. The reconstrucrion 
required only 6 bytes and represented a compression to 2 
percent of the original. In this case stage two compression 
was not applied because it increased the number of bytes to 14. 

Figs. 14A-E illustrates the trade between compression 
ratio and the fidelity of the reconstruction. Fig. 14A is the 
original of a printed address which required a total of 1808 
bytes. Fig. 14B shows the first level of compression which 
employed a compression ratio similar to that used in Figs, lie 
and 12C. After both stage one and stage two compression, this 
reconstruction is 13 percent of the original. Figs. 14C-E show 
gradually increasing levels of- compression. The reconstruction 
shown in Fig. 14C, requiring only 9 percent of the original 
appears quite similar to the original. Figs. 14D-E, although 
quite readable, show definite differences from the original 
The reconstructions in Figs. 14D and 14E require 8 percent and 
7 percent of the- original storage space, respectively. 
Although some applications require high fidelity, such as 
signature verification, many applications could tolerate the 
lower fidelity accompanying the high compression ratios shown 
m Figs. 14D-E. 

Fig. 15 shows an embodiment of the apparatus to 
accomplish the method described in this invention. in the 
preferred embodiment, the computer 1550 is an IBM/PC/AT 
compatible with an IBM compatible monitor 1540 and keyboard 
1510. computer 1550 employs the methods disclosed by this 
invention to compress the electronic ink. Handwritten 
information is input either directly via the hardwired 
digitizer tablet 1520 or via a portable digitizer tablet 1530 
Portable digitizer tablets are in common use in many 
industries, including the package delivery industry where an 
electronic signature is obtained upon delivery of a package. 
Depending upon the computing power of tablet 1530, the 
signature can either be immediately compressed using this 
invention's compression methods or temporarily stored for 
subsequent downloading to computer 1550. Once the handwritten 
aata has been compressed, it can either be stored in computer 
looO or in separate memory 1560. Monitor 1540 allows the 
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* original handwritten data to be compared with the compressed 
data. 

The invention has now been explained with reference to 
specific embodiments, other embodiments will be apparent to 
those of ordinary skill in this art in light of this 
disclosure. It is therefore not intended that this invention 
be limited except as indicated by the appended claims. 
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APPENDIX 

The computer program modules used to compress the 
5 elecrronic data are given in this appendix. Significant 

program modules and their specific functions are as follows: 

Stage one compression: 

calc_sav.c Function to determine which points in the 
10 signature need to be saved, based on various 

criteria. 

curvature. c Function to calculate the local curvature at a 
point along the contour. 

Stage rwo compression: 

cmpsig2.c Function to compress an "array of stage 1 data 

with characteristics of the signature 
verification template data. 
pack_half bytes . c 

Function to pack the value to be compressed into 
half bytes. 

set_4_bits Function to pack values into the upper or 

lower 4 bits of a byte. 
gst_4_bits Function to unpack values from the upper or 

lower 4 bits of a byte. 
uncmpsig2.c Function to uncompress an array of data 

compressed by the function cmpsig2.c. 
unpack_half bytes. c 

Function to unpack and reconstruct the value 
from half bytes. 
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35 

/* CALC^SAV.C 

★ ^. , ^ 

^ Function zo determine which points in the signature need to 
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* be saved, based on various criteria. The results are 

* returned in array save_ptsc. 

* For example, if the ith point of x[],yn should be saved 

* then save_ptsC[i] = i, else save _pts[i] = O; 

* cur.comp returns: i (i.e. yes or TRUE) for successful 

scaling 

0 (i.e. NO or FALSE) if there is an 
error 



input arguments: curvature.threshold, nsamps, 



segpnt, x, y 



* output arguments: save_pts 



numsegs, 



^include <stdio.h> 
#include <cic.h> 

BOOL calc_sav(curvature_threshold, nsamps, numsegs, segpnt, x, 
y, save_ptsC) ' ' 



COUNT curvature_threshold; 
COUNT nsanps; 

COUNT nuir.segs; 

COUNT segpnt [ ] ; 

COUNT X [ ] ; 

COUNT y [ ] ; 

TINY save_ptscr ] ; 



/* threshold for curvature*/ 



COUNT i; 
COUNT k; 



/* compression 
/* number of samples 
/* (length) in input 
/* arrays x and y 
/* number of individual 
/* segments (strokes) 
/* pointers to beg of 
/* individual strokes 
/* input array of x 
/* values 

/* input array of y 
/* values 

/* contains information 
/* about what points to 
/* keep (returned to 
/* calling program) 
/* loop counter 
/* loop counter 



*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 



WO.94/038S3 

COUNT valuel; 
COUNT value2; 
COUNT checks; 
COUNT lower_liinit ; 

5 

COUNT upper^limit ; 
BOOL curvature ( ) ; 
0 COUNT si2e_angle; 

BOOL clockviseB; 
BOOL calc_dir; 

5 

COUNT dir_angle; 
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/* 


temporary variable 


*/ 


/* 


temporary variable 


*/ 


/* 


temporary variable 


*/ 


/* 


defines starting pt 


*/ 


/* 


for curvat compress 


*/ 


/* 


defines ending pt for 


*/ 


/* 


curvat compress 


*/ 


/* 


function to calculate 


*/ 


/* 


local curvature 


*/ 


/* 


calculated curvature 


*/ 


/* 


at a point 


*/ 


/* 


YES if clockwise, NO 


*/ 


/* 


otherwise 


*/ 


/* 


if YES calculate 


*/ 


/* 


dir_angleP 


*/ 


/* 


direction of curvature 


*/ 


/* 


bisector 


*/ 



* Initialize the save_ptsc array. 
0 * . 

*/ 

for (i = 0; i < nsamps ; -+i) 

{ 

save_ptsC[i] = 1; 

=> ) 
/* 

■«■ 

" Perform curvature based compression (determines which points 
* to keep of the stroke and saves the information in array 
0 * save_ptsC to return to the calling program) . 

•K — — — — — — — — — 

for (i = 0; i < numsegs; +-t-i) 

I 

5 lower_liinit = segpnt[i]; 

upper_liinit = segpnt[i + i]; 
while (lower^limit < upper_liinit ) 



10 



ID 



20 



25 



30 
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checkB = NO; 

for (k = lower^liinit + 2; k < upper^limit; k += 2) 

( 

curvature (X, y, (k + lower^liitiit) / 2, 
nsamps, (k - lower^limit) / 2, 
&si2e_angle, &clockwiseB, &calc_dir, 
&dir_angle) ; 
if (size_angle > curvature_threshold) 
{ 

CheckB = YES; 
break; 

} 

) 

if (CheckB) 

i 

value2 = (k + lower^limit ) / 2; 
valuel = value2; 

) 

else 

{ 

valuel = segpnt[i +1] - i; 

value2 = (valuel + lower^limit) / 2; 

) 

for (k = lower_limit 1; k < valuel; — rk) 

/ 

if (k != value2) 

save_ptsc[k] = O; 

) 

if (CheckB) 

lower_limit = value2 ; 

else 

lower_liinit = upper_limit; 

) 



/ 



* 



T to calling program with successful completion fl 
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* ^ 

*/ 

return ( YES ) ; 

/* 

* 

*/ 

J /* end of function */ 

/* CALC_SAV */ 

/* CURVATURE. C 
* ' 

* Calculates the curvature, direction of curvature (angle) 

* bisector, and whether the change in angle is clockwise or 

* counter_clockwise at a point defined by function argument 

* "index" into the arrays of coordinates points x(], y[], 

* representing the stroke contour. The quantities returned 

* are in tenr.s of the angle range 0-360 degrees. 
* 

* If ERROR is detected the BOOL function value is YES (=1) , 

* while for a normal return the BOOL function value is NO 

* (=0). 

★ 

* CAUTION: when rhe curvature is 0 (i.e., calculated 

- si2e_angleP = 0 or close to 0), the calculation of the 

* direction of the curvature bisector is, of course, 

* inherently ambiguous and unreliable. 
•k . ^ 

* / 

#include <stdio.h> 
^include <cic.h> 



3C0L curvature (X, y, index, ma 
clockviseBP, calc_di 
COUNT X [ J ; 

COUNT y [ 1 ; 

COUNT index; 



xpts, Ighseg, size_angleP, 

.rB, dir_angleP) 

/* X coordinate values of */ 
/* character */ 
/* y coordinate values of */ 
/* character */ 
/* point in x, y to calc. */ 
/* curvature */ 
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/* 



COUNT maxpts; 



COUNT Ighseg; 



COUNT *size_anglep; 



BOOL *clockwiseBP; 



BOOL calc_dirB; 



COUNT *dir_angleP; 

{ 

COUNT plus_point; 
COUNT rainus_point; 
COUNT plus_angle; 

COUNT minus_angle; 

COUNT dif_angle; 

COUNT cic_atan2 0 ; 



*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 



/* maximum length of 
/* arrays x,y 
/* length of segment for 
/* curvatur calc 
/* calculated curvature 
/* at point index 
/* YES if clockwise, NO 
/* otherwise 
/* if YES calculate 
/* dir_angleP, otherwise 
/* otherwise return 
/* dir_angle = o 
/* direction of curvature */ 
/* bisector */ 

/* stores "index + Ighseg"*/ 
/* stores "index - Ighseg"*/ 
/* stores calc direction */ 
/* of upper seg. 
/* stores calc direction 
/* of lower seg. 
/* stores angle angle 
/* difference 
/* function to calculate 
/* angle in the 
/* range -1800 to 1800 
/* (degrees x 10) 



*/ 

*/ 

*/ 

*/ 
* / 

*/ 
*/ 
V 
*/ 



* Check for errors in specifying the range for 

* computation. 



curvature 



minus_point = index - Ighseg; 
plus_point = index ^ Ighseg; 
if (ininus_point < o 1 [ plus_point >= maxpts) 



*si2e_angleP = o; 
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*clockwiseBP = YES; 
*dir_angleP = 0; 
return (YES) ; 



/* ERROR - out of bounds, */ 
/* cannot calc. curvature */ 



* Calculation of the local curvature (i.e., angle change). 

* Note: the function cic_atan2 returns angle values in the 

* range -1800 to 1800, so it is necessary to divide by lo to 

* get the desired curvatures. 



* / 



plus_angle = cic_atan2 (x(plus_point] - x[index], 

y[plus_point] - y[index]); 
minus_angle = cic_atan2 (x ( index] - x[ininus_point] , 

y( index] - y [niinus_point] ) ; 
dif_angle = plus_angle - ininus_angle; 
if (dif_angle > 1800) 

dif_angle -= 3 600; 
else if (dif_angle < -1800) 

dif_angle += 3 600; 
if (dif_angle < 0) 

( 

*size_angleP = (-dif_angle 5) / lO; 
*clockwiseBP = YES; /* by convention 

/* clockwisfi 
' /* angle changes are 

/* negative 

else 



V 
*/ 
*/ 
*/ 



*si2e_angleP = (dif_angle + 5) / lo; 
*clockwiseBP = NO; /* by convention 

/* counter-clockwise 
' /* angle changes are 

/* positive 



*/ 
*/ 
*/ 
*/ 



10 



15 



0 
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* Calculate the direction of the bisector of the local 

* curvature (i.-e., angle change) computed above in the range 

* 0-360 degrees. Note: the cic_atan2 returns angle values in 

* the range -1800 to 1800, so an adjustment is made in the 

* code below so that the resulting calculated bisector value 

* is in the range 0-3 60. 
* — 

*/ 

if (calc_dirB) 

{ 

if (*clockwiseBP == YES) 

*dir_angleP = plus_angle - (900 + (dif_angle - i) 
/ 2) ; 

else 

*dir_angleP = plus_angle + (900 - (dif_angle + i) 
/ 2) ; 

if (*dir_angleP < O) 

*dir_angleP += 3 600; 
*dir_angleP = (*dir_angleP +5) / lO; 
if (*dir_angleP 360) 

*dir_angleP -= 3 60; 

) 

/* normal return - no */ 
/* errors detected */ 



*/ 
/* 

0 CMPSIG2.C 



/* end of function 



/ 



/* CURVATURE. C */ 



* Function to compress an array of data with charcteristic 

* the signature verification template data. The companion 

* program to uncompress the data is UNCMPSIG2.C. NOTE: 
compression/uncompression is exact (not lossy) and 

* restores the original data exactly. 



* 



cmpsig2 returns: 1 (i.e. yes or TRUE) for successful 
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* compression 

* 0 (i.e. NO or FALSE) if there is an 

error 



* input arguments: 

* numin - number of elements (length) of data array to 

* be compressed 
array - data array to be compressed 

maxout - maximum allowed size of the compressed array 
comp_arrayUC[ ] (to prevent overwriting 
memory) 



* 
* 

10 * 



15 * 
* 

20 * 



* ourput arguments: 

* num_comp_bytesP - ptr to # elements (each a byte) of 

the compressed data stored in 
comp_arrayUC[ ] . 

comp_arrayUC - unsigned byte array containing the 
compressed representation of the input 
data in array [ ] ; 



#include <stdio.h> 
=include <cic.h> 



25 BOOL cmpsig2 (numin, array, : 

comp_arrayUC) 
COUNT numin; 

COUNT array [ ] ; 

30 

COUNT maxout; 
COUNT *num_comp_bytesP ; 
35 TBITS comp_arrayUCr ] ; 

{ 

COUNT i; 



, num_comp_bytesP , 



/* number of points in input */ 

/* array ★/ 

/* array containing input */ 

/* data to be compressed */ 

/* max number of pts for */ 

/* comp^arrayUC ★/ 

/* actual #bytes used in */ 

/* comp_arrayUC */ 

output array of compressed */ 

/* data */ 

/* loop index */ 
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COUNT num_sign_change; /* counts number of sign */ 

/* changes */ 

COUNT value; /* for general use temporary */ 

/* storage */ 

COUNT byte_part; /* 1 if upper 4 bits of byte, */ 

/* and 2 if lower 4 bits of */ 
/* byte */ 

TBITS *comp_arrayUCP; /* ptr to comp^arrayUC of sig.*/ 

/* data 

TBITS byte_maskUC; /* byte used to temporarily */ 

/* hold compressed data */ 

COUNT last^index; /* saves last index of direct.*/ 

/* change 

COUNT current^sign; /* used in sign change */ 

/* determination ★/ 

COUNT last_sign; /* used in sign change */ . 

/* determination */ 

BOOL set_4_bits() ; /* function to pack values */ 

/* into the upper or lower */ 
/* 4 bits of a byte */ 

BOOL get_4_bits() ; /* function to unpack values */ 

/* from the upper or lower */ 
/* 4 bits of a byte */ 

BOOL pack_halfbytes() ; /* implements main packing */ 

/* logic 

/* 

★ ^ - ^ 

* Initializations for the compression process. 
* 

*/ 

*nuiii_coinp_bytesP = O; 
*corop_arrayUCP = conp_arrayUC ; 

/* 

* First, save the number of items in the array to be 

* compressed. This value is stored as two bytes. The first 

* byte is the upper 8 bits of the 16 bit COUNT value for 

* numin, and the second byte is the lower 8 bits. 
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★ ^ ^ ^ ^ 

*/ 

*nuin_coinp_bytesP += 1; 

*coinp_arrayUCP++ = (TBITS) ((numin >> 8) & Oxf f ) ; 
5 *nuin_coinp_bytesP += 1; 

*coinp_arrayUCP++ = (TBITS) (numin & Oxff ) ; 

/* 

* Find the number of sign changes (i.e., when successive 
0 * difference values have different signs) and save to 

* coinp_arrayUC. Note that the maxiinuKi number of sign changes 

* allowed is 255. The function immediately exits and returns 

* an error flag if > 255 sign changes are detected. 

5 */ 

num_sign_change = 0; 

last_sign = array[l] - array[0]; 

for (i = 2; i < numin; ++i) 

i 

0 current_sign = array[i] - array [i - 1]; 

if ( (current_sign < o && last_sign >= 0) 1] 
(current_sign >= o last^sign < 0) ) 

{ 

last_sign = current_sign ; 
5 •^+num_sign_change ; 

) 

) 

if (num_sign_change < 2 56) 

{ 

0 *num_comp_bytesP 1; 

*comp_arrayUCP-H+ = (TBITS) num_sign_change ; 

) 

else 

{ 

5 fprintf (stderr, "\n\nCMPSIG — num_sign_change > 

256\n\n") ; 
return (NO) ; 
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r. <» 

/* 

* — - — 

* Initializations for compression 

* . 

*/ 

byte_part = 1; 
byte_maskUC = 0; 

/* 

» 

* Pack the sign of the first value of the input data (i.e., of 

* array(0]) 

* 

V 

if {array[01 >= 0) 
value = 1; 

else 

value = 0; 

pack_half bytes (value, &byte_part, &byte_inaskUC, 
&comp_arrayUCF, nuin_coTnp_bytesP) ; 



* Pack the sign of the first difference value (i.e., of 

* array [11 - array [0]) 



*/ 

last_sign = array [1] - array [ 0]; 
if (last_sign < 0) 
value = 0; 

else 

value = 1; 

pack_halfbytes( value, &byte_part, &byte_maskUC, 
&coinp_arrayUCP, nuin_coHip_bytesP) ; 

/* 



* Pack the sign change positions (where successive difference 

* values have different sign) . 
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last_index = 0; 

for (i = 2; i < numin; ++i) 

{ 

current^sign = array [i] - array [i - i]; 
if ( (current_sign < o && last^sign >= O) |j 
(current^sign >= o last^sign < o) ) 

{ 

last_sign = current_sign; 
if ((i - last_index) > 255) 

{ 

fprintf (stderr, "\n\nCMPSIG ~ sign index 

change > 255\n\n") ; 
return (NO) ; 

) 

pack_halfbytes{i last^index, &byte_part, 
&byte_maskUC, &comp_arrayUCP, 
nuii_coinp_bytesP) ; 

last_index = i; 

) 

^ /* end of loop to pack sign */ 

/* change position values */ 



- Pack differences values into bytes and save to coinp_arrayUC 

*/ 



for (i = 0; i < numin; 
{ 

if (i == 0) 

value = abs (array [0] ) ; 

else 

value = abs (array [i] - array [i - i]); 
pack_half bytes (value, &byte_part, &byte_naskUC, 

&coinp_arrayUCP, nuni_conip_bytesP,). ; 
if (*num_comp_bytesP >= maxout) 

{ 

fprintf (std-rr,"\n\nCMPSIG — conip_arrayUC too 
large\n\n") ; 
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5 /* 



return (NO) ; 

) 

/* end of loop to pack ★/ 
/* difference values */ 



* Finally, if the data ends on an odd number of half bytes, 

* must save one final byte because the last signature data 

* value is in the first half byte of byte^maskUC and hasn't 
10 * yet been saved to comp^arrayUC. 

* ^ ^ 

V 

if (byte_part == 2) 



set_4_bits((UTINY) 15, 2, &byte_inaskUC) ; 

*nuin_coinp_bytesP += 1; 

if ( *nuin_coinp_bytesP > maxout) 

{ 

fprintf (stderr, "\n\nCMPSIG comp^arrayUC too 
20 large\n\n") ; 

return (NO) ; 

) 

*coTnp__arrayUCP++ = byte_maskUC; 

) 

25 /* ' 

* 

"^is point is reached the compression process he 



las 

* proceeded normally and no e.-rors have been detected, so 

* return with a "successful completion" status flag. 



30 *• 

rerurn(YES) ; 
★ . ^ 

V 

35 1 

' /* end of function CMPSIG2 */ 

/* ?ACK_HALFBYTES , C 

•K — —— — — — — — 

* Function to pack a COUNT value into half bytes. 
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* 
* 

* 

*/ 



pack_half bytes returns: 1 (i.e. YES or TRUE) for 

successful packing 
0 (i.e. NO or FALSE) if 
there is an error 



BOOL pack_half bytes (value, byte_partP, byte_iaaskUCP, 
coinp_arrayUCPP, nuni_comp_bytesP) 

COUNT value; /* value to be compressed */ 

COUNT *byte_partP; /* i if upper 4 bits of byte, */ 

/* and 2 if lower 4 bits of */ 

/* byte */ 

TBITS *byte_inaskUCP; /* packed byte representing */ 

/* different values */ 
TBITS **coinp_arrayUCPP; /* ptr to ptr to comp.arrayUC */ 

/* element */ 

COUNT *nuin_coinp_bytesP; /* actual kbytes used in */ 

/* conp_arrayUC */ 

UTINY tempUC; /* temporary value */ 

UTINY continuation_labelUC;/* used to save continuation*/ 

/* label */ 

BOOL get_4_bits() ; /* function to unpack values */ 

/* from the upper or lower 4 */ 

/* bits of a byte */ 

BOOL set_4_bits() ; /* function to pack values */ 

/* into the upper or lower 4 */ 

/* bits of a byte */ 

if (value < 14) 

{ 



it . 



* Pack 4 bit value into byte. When the absolute difference 

* value is < 14, then it can be packed into a half_bvte (4 

* bits) . 



*/ 
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set . bits,, CTINV, value, .byte_partP, byte..asWcP, 
If (*byte_partP ==2) 

{ 

*byte_partP = i- 

*nuin_conip_bytesP += i ; 
**coinp_arrayUCPP = *byte_maskUCP; 
++(*coinp_arrayUCPP) ; 

) 

else 

*byte_partP = 2; 

) 

else 



. IITJ. =°"^~"^=n label into byte. When the absolute 
difference value is <. 30 (and, of course, >= u since the 
previous block above deals with values < 14, ILl 

* value 14 ir,^-« aj-ues < 14), then pack a 

. K ^""^ 'he re,,ainder 

into a subsequent half byte. If the absolute difference 
value .s >. 30. pacK 15 into the first half byt^ and acK 
th re.a„der into another byte. Zf the absolute J fe ence 

value IS > 14 + 255 = -x^Q -^^J-eirence 

* ' ^" has occurred and 

the function immediately exits anrf burred and 

* condition. ^ ^^^^^"^ 



* 



/ 



if (value > 269) 

{ 

fprintf(stderr, "\n\nCMPSlG - abs dif value out of 

fprintf(stderr, "range (= %d)\n\n", value)- 
return (NO); 

^ ' >' /* return on ERROR detection */ 

else if (value < 30) 

continuation_labelUC = (UTINY) 14; 

else 

continuation_labelUC = (UTINY) 15; 
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/* 

★ . ^ «. « 

* First, set half byte for continuation label. 



0 



7 








set_4_bits (continuation^ 


^labelUC, *byte_partP, 




byte^maskUCP) ; 






if (*byte_partP == 2) 

( 






*byte_partP 


1; 




*nuin_comp_bytesP 


1; 




^ wAii^^Q i. i a y u — 


*Dy te_iBasKUCP ; 




(*coinp_arrayUCPP) ; 

) 






else 






*byte_partP = 2 ; 






if (value < 30) 

( 





* Pack a value corresponding to (value - 14) into a second 

* half byte (i.e., a value from 14 to 29 will be converted to 

* a value f rcn 0 to 15 for packing) . 
* 



set_4^bits( (UTINY) (value - 14), *byte^partP, 

byte_maskUCP) ; 
if (*byte_partP == 2) 

{ 

*byte_partP = i; 

*nuin_coinp__bytesP -;-= l; 
**comp_arrayUCPP = *byte_inaskUCP ; 
( *coinp_arrayUCPP) ; 

) 

else 

*byte_partP = 2 ; 
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/* 



* If the original absolute difference value is >= 3o'"thir' 

* must pack two more half bytes in additional to the 

* continuation half byte. 

V 

if (*byte_partP == i) 

' { 

*byte_inaskUCP = (tbITS) (value - 14); 

*nuin_coinp_bytesP += i; 
**coinp_arrayUCPP = *byte_inaskUCP; 
++ ( *conip_arrayUCPP) ; 

) 

else 

{ 

get_4_bits(&tempUC, l, (tbiTS) (value - 14)). 
set^4_bits(teinpUC, 2, byte_niaskUCP) ; 
*num_comp_bytesP += i; 
**coinp_arrayUCPP = *byte_inaskUCP; 
++ ( *coinp_arrayUCPP) ; 

get_4_bits(&teinpuc, 2, (TBITS) (value - 14)). 
set_4_bits((UTINY) tempUC, i, byte_naskUCP) ; ' 

) 

) 

return (YES) • 

^^^"e unpacked successfully*/ 

/* end of pack halfbvtes */ 
/* raCMPSIG2.C -mxroyres */ 

* Function to uncompress an array of data comprlslld'by 

.unc.on CMPSIG2. NOTE: co.pression/unco.pression is exact 
(not lossy) and restores the original data exactly 



uncr,psig2 returns: i (i.e. yes or TRUE) for 

successful uncompress 
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* 0 (i.e. NO or FALSE) if 

* there is an error 
* 

* input arguments: 

5 * comp^arrayUC - unsigned byte array containing the 

* compressed representation of the data 

* (as created by CMPSIG2.C) 
★ 

* output arguments: 

10 * numoutP - ptr to number of elements (length) of the 

* uncompressed data 

* array - data array created by unpacking and 

reconstructing values taken from the 
compressed array comp_arrayUC [ ] . 



★ 

15 * 



#include <stdio.h> 
^include <cic.h> 



2 0 300L uncnpsig2 (numoutP, ai 

COUNT *numoutP; 

COUNT array [ j ; 

25 

TBITS comp_arrayUC[ ] ; 

( 

3 0 COUNT i; 

COUNT kount^bytes; 

TBITS *comp_arrayUCP; 
3 5 COUNT num_sign_change; 

COUNT value; 



, comp_arrayUC) 



/* 


number of points in output 


*/ 


/* 


array 




*/ 


/* 


array for output signat 


*/ 


/* 


data that has reconst 


from 


*/ 


/* 


coinp_arrayUC by reversing 


V 


/* 


the compression proc. 




V 


/* 


compressed input array 


of 


*/ 


/* 


sig data 




*/ 


/* 


loop index 




*/ 


/* 


kounts bytes used in 




*/ 


/* 


comp_array 




*/ 


/* 


ptr to comp_arrayUC of 


sig 


*/ 


/* 


data 




*/ 


/* 


counts number of sign 




*/ 


/* 


changes 




*/ 


/* 


for general use temporary 


*/ 


/* 


storage 




*/ 



10 



15 



20 



25 



30 



35 
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COUNT byte_part; /* i if ^ppe^ , bits of byte, */ 

/* and 2 if lower 4 bits of */ 
/* byte 

TBITS byte.maskUC; /* packed byte representing */ 

/* different values */ 

COUNT sign_change[256]; /* array of sign change */ 

/* indices 

COUNT *sign_changeP; /* pointer to sign change */ 

/* index array */ 

COUNT sign_first_dif; /* saves sign of first */ 

/* difference */ 

COUNT current_sign; /* current value of sign */ 

BOOL unpack_halfbytes(); /* implements main unpacking */ 

/* logic */ 

/* 

■k 

* Initializations for compressed array 

-k 

*/ 

kount_bytes = O; 
comp_arrayUCP = comp_arrayUC; 

/* 

* 

* First, get the number of items in the array that was 

* compressed. This value is stored as two bytes. The first 

* byte is the upper 8 bits of the 16 bit COUNT value for 

* numout, and the second byte is the lower 8 bits. 

*/ 

*numoutP . = 0; 

*numoutP = *comp_arrayUCP++; 

*numoutP <<= 8; 

*numoutP = *numoutP I ( *comp_arrayUCP+^) ; 

kount_bytes += 2 ; 
/* 

★ 



* Get the number of sign changes 

X — — — — _ — — 
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*/ 



num_sign_change = (COUNT) (*coinp_arrayUCP++) ; 
++kount_bytes ; 

/* 

■k 

* Check for error condition. If the value rlld'in^or'thr 
number of sign changes is >= 256, then an error has 

* occurred. 



*/ 



if (num.sign.change < o | | nuin_sign_change >= 256) 

fprintf(stderr, "XnVnUNCMPSIG - number of sign 

changes out " ) ; 
fprintf(stderr, -of range (= %d) is >= 256\n\n", 

nuin_sign_change) ; 
return(NO) ; ^^^^ detection */ 

/* 

* 

* Initializations for unpacking. 

byte_parr = i; 

/* 

^ 

* unpac. the sign of the first value^rthrori^Inai" 

* uncompressed data. 
* 

V ■ 

unpack.halfbytes (.value, &byte_part, .byte_maskUC, 

S.comp_arrayUCP, &kount_bytes) ; 
if (value == o) 

current_sign = -i; 
else if (value == i) 

current_sign = i; 
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/* 

* 

* Error condition since, value != 0 (for roinus sign) and also 

* != 1 (for plus sign). 

5 *— 

*/ 

fprintf (stderr, "\n\nUNCMPSIG — initial sign value ••) ; 

fprintf (stderr, "is inconsistent (= %d)\n\n", value); 
return (NO) ; /* return on ERROR detection */ 

0 ) 

/* 

★ . ^ 

* Unpack the sign of the first difference value of the 

* original, uncompressed data. 

5 * 

*/ 

unpack^half bytes (&value, &byte_part, &byte_inaskUC, 
fitcomp^arrayacp, 
&kount_bytes) ; 
) if (value == O) 

sign_f irst_dif = -1; 
else if (value == 1) 

sign_f irst_dif = l; 

else 

^ { 
/* 

★ ^ 

* Error condition since value != 0 (for minus sign) and also 

* i= 1 (for plus sign) . 

\ *- — 

V 

fprintf (stderr, "\n\nUNCMPSIG — initial sign value ") ; 
fprintf (stderr, "is inconsistent (= %d)\n\n", value); 
return (NO); /* return on ERROR detection */ 

) 

/* 

* 



* Unpack sign change positions (where successive diff 



erence 
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* values have different sign) 
* 

*/ 



for (i = 0; i < nuin_sign_change ; ++i) 

5 { 

unpack_halfbytes(&value, &byte_part, &byte_inaskUC, 
&comp_arrayUCP, &kount_bytes) ; 

if (i == 0) 

sign_change[0] = value; 

10 else 

sign_change[i] = sign_change[ i - 1] + value; 
J /* end of loop to unpack sign */ 

/♦change position values */ 



20 



25 



30 



/* 



* Unpack differences from bytes and read from disk. 



*/ 



sign_changeP = sign_change ;/* pointer to sign change */ 

/♦index array */ 
for (i = 0; i < *numoutP; ++i) 

{■ 

unpack_half bytes (Svalue, &byte_part, &byte_inaskUC, 

&comp_arrayUCP, &kount_bytes) ; 
if (i > 1 && nura_sign_change > 0 && i « *sign_changeP) 

{ 

current_sign = -current_sign; 
— num_sign_change ; 
++sign_changeP ; 

) 

if (i == 0) 

array [ii = current_sign * value; 
else if (i == i) 

{ 

sign_cha ? = sign_change; 

arrayfi; = array [i - 1] + sign_f irst_dif * value; 

current_. gn = sign_f irst_dif ; 
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else 

array [i] = array [i - 1] + current_sign * value; 
) /* end of loop to unpack and */ 

/* reconstruct data */ 

5 /* 



* Check for inconsistency in specified size (in bytes) of the 

* input compressed array {coinp_array) versus the actual number 

* of bytes used in the uncompress process to create the output 
10 * uncompressed array. 

★ «^ 

if (kount_bytes 1= num_comp_bytes) 

{ 

fprintf (stderr, "\n\nUNCMPSIG ~ kount^bytes (= %d) 

i= " , kount_bytes) ; 
fprintf (stderr, "num_comp_bytes (= %d)\n\n", 
num_comp_bytes) ; 

return (NO) ; 

) 



* If this point is reached the uncompression process has 

* proceeded normally and no errors have been detected, so 

* return with a "successful completion" status flag. 



V 




return (YES) ; 








30 ) 


/* end of function UNCMPSIG2 */ 


/ ^ UNPACK_HALFBYTES , C 




* Function to unpack 


and reconstruct a COUNT value from half 


* bytes. 




35 * 




* pack_halfbyres 


returns: 1 (i.e. YES or TRUE) for 


* 


successful unpacking 


* 


0 (i.e. NO or FALSE) if 



20 */ 
/* 
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there is an error 



*/ 

BOOL unpack_hal f bytes (valueP, byte_partP, byte^maskUCP, 
comp_arrayUCPP, kount_bytesP) 



COUNT *valueP; 

COUNT *byte_partP; 

TBITS *byte_iriaskUCP; 
TBITS * *conip_arrayUCPP ; 
COUNT *kount_bytesP; 



/* contains uncompressed 
/* value to be returned to 
/* main program 
/* 1 if upper 4 bits of byte,*/ 
/* and 2 if lower 4 bits of */ 
/* byte 

/* packed byte representing 
/* different values 
/* ptr to comp_arrayUC of sig 
/* data 

/* kounts bytes used in 
/* comp^array 



*/ 
*/ 
*/ 



*/ 
*/ 
*/ 



*/ 
*/ 
*/ 
*/ 



UTINY tempUC; 
2 0 UTINY valueUC; 

BOOL get_4_bits() ; 



BOOL set_4_bits() 



/* temporary value 

/* temporary value 

/* function to unpack values 

/* from the upper or lower 4 

/* bits of a byte 

/* function to pack values 

/* into the upper or lower 4 

/* bits of a byte 



it. 



*/ 
*/ 
*/ 
*/ 
*/ 
V 
*/ 
*/ 



* unpack 4 bits. When the value is < 14, it corresponds to a 

* legitimate difference value, and hence the next array value 

* can be calculated by adding the difference (with the correct 

* sign) to the preceding array value. If the value == 14, ■ the 

* absolute difference value is >= 14 and <= 29, hence another 

* half byte must be read in to compose the absolute 

* difference. if the value == 15, the absolute difference 

* value is >= 30 and < 269, hence another full byte (i.e., 2 

* half bytes) must be read in to compose the absolute 

* difference. 
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— — — - — — - — — — — — — — « 

*/ 

if ( *byte_partP == 1) 

{ 

++ ( *kount_bytesP) ; 

*byte_raaskUCP = **comp_arrayUCPP; 

++ ( *coinp_arrayUCPP) ; 

) 

get_4_bits(&valueUC, *byte_partP, *byte_niaskUCP) ; 
*byte_partP = 3 - *byte_partP; 
if (valueuc < 14) 

*valueP = (COUNT) valueUC; 
else if (valueUC == 14) 

{ 

if (*byte_partP == l) 

( 

*byte_maskUCP = **comp_arrayUCPP; 
++(*coinp_arrayUCPP) ; 
++(*kount_bytesP) ; 
) 

get_4_bits(& valueuc, *byte_partP, *byte_maskUCP) ; 
*valueP = (COUNT) (valueUC +14); 
*byte_partP = 3 - *byte_partP; 

) 

else if (valueUC == 15) 

( 

if (*byte_partP == i) 

{ 

. *byte_maskUCP = **conip_arrayUCPP; 
++ ( ''Coinp_arrayUCPP) ; 
++(*kount_bytesP) ; 
) 

get_4_birs(&t6:npuc, *byte_partP, *byte_maskUCP) ; 
set_4_bits(tempUC, 1, &valueUC) ; 
*byte_partP = 3 - *byte_partP; 
if (*byx:e_partP == i) 
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*byte_maskUCP = **conp_arrayUCPP; 
++(*coinp_arrayUCPP) ; 
++ ( *kount_bytesP) ; 

) 

get_4_bits(&teinpUC, *byte nartP *k„i. 

«yte_partp, *byte_inaskUCP) ; 

set,4_bits(teinpuc, 2, ivalueUC) ; 
*valueP = (COUNT) (valueUC + 14). 
*byte_partP = 3 - ^bytejartpj 

else 

( 

*byte_«,askUCP = **comp_arrayUCPP; 
++(*coinp_arrayuCPP); +.-(*kount bytesP) • 

return (YES) ; /* 
/ ""^^"s unpacked successfully */ 

*/ 

\ /* end of function UNPACK_HALFBYTES */ 

/* GET_4_BITS.C 
*-----—-———__ 

* Function that returns the UTINY val„o 

* eithPr ^ho ^ corresponding to 
ej-tner the upper 4 hii-c •;„ t-i. • 

* bits Of ^"P"'' '=1"=^' "'^ -hs lower « 
firs. Of course this be a value in 

* The input arguemeht "pcition.. a^Z, ! ° " 

* or lower 4 bits is taLn: '"""""^^ ""ether the upper 

position = 1 ^ 

upper four bits (or upper 

^ "nybble") 

position = 2 1^.,^ ^ 

lower four bits (or upper 

"nybble") 



??include <stdio.h> 
=include <cic,h> 

JXm .valueUCP; ^^^^ ^.^^ 

/* value 



V 

*/ 
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COUNT position; , ' ' 

/* 1 specifies the upper 4 */ 

/* bits, and ^. 
TBITS byte maskUC; /* ^^,^„ , ^ 

/* feature mask for template */ 

/* selection 2 specifies the */ 

/* lower 4 bits 
{ */ 

/* 

* 

co'ndi'tiln'fo'und 

*/ 

if (position 1= 1 i& position != 2) 
return (NO) ; 

/* 

★ ^ 

*/ 

if (position == i) 

byte^maskUC >>= 4; 
*valueUCP = (UTINY) (byte_maskUC & Oxf ) ; 

return (YES) ; ^* "^^^^ V 

* ^ 

V 

*/ ' /* ^"'^ °* function GET_4_BITS.CV 

* SET_4_BITS.C 
•* 

* Function that sets either the uppIr^^bitrinTh""' 

* bvte or 1 input 

^ specified by the input argument "value". 
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* "position =1 — upper four bits (or upper 

* "nybble") 

* position -2 — lower four bits (or upper 

* "nybble") 
5 * 

* NOTE: function returns a boolean YES for normal operation 

* and NO if an error is detected. 

» 

*/ 

0 include <stdio.h> 
iinclude <cic.h> 

BOOL set_4_bits (valueUC, position, byte_inaskUCP) 

UTINY valueUC; /* set the four bits to this */ 

/* value */ 
COUNT position; /* specifies the position of */ 

/* the two bits */ 
TBITS *byte_naskUCP; /* feature mask for template */ 

/* selection */ 

{ 

) /* 

•k — 

Return original mask if error condition found 

ir — — — — — — — — — — 

V 

> if (valueUC > (UTINY) 15) 

return (NO) ; 
if (position != l && position != 2) 
return (NO) ; 

* Set bits as specified by the input arguments 

— — — — — — — — — 

*/ 

if (position == i) 

{ 

*byte_maskUCP &= "(OxfO); 

/* zero upper 4 bits inbyte_inaskUCP*/ 
*byte_maskUCP |= (valueUC << 4); 
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/* set upper 4 bits to new value */ 
else 

{ 

5 *byte_maskUCP &= •(Oxf); 

/* zero lower 4 bits in byte_riaskUCP*/ 
*byte_maskUCP j= valueUC; 

/* set upper 4 bits to new value */ 

10 return(YES); 

/* 

^ /* end of function SET 4 BITS * / 

/* - - . / 



15 * This file contains standard type definitions, defined 
* constants, macros, and data structures. 
* 

*/ 

#define FAST register 
20 #define GLOBAL 

=define IMPORT extern 

=define THISFL extern 

#define INTERN static 

#define LOCAL static 

25 

typedef void VOID; 

typedef char TBOOL, TINY, TEXT; 

typedef short COUNT, SHORT; 

typedef int int, ARGINT, METACH; 

3 0 typedef long LONG; 

typedef unsigned char UTINY, TBITS, UTEXT; 

typedef unsigned short BITS, UCQUNT, USHORT, SIZETYPE; 

typedef unsigned int BYTES, UINT; 

3 5 typedef unsigned long ULONG, LBITS ; 



typedef float FLOAT; 
typedef double DOUBLE; 
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typedef USHORT ERCTYPE; 

typedef enum 

{ 

5 NO = 0, 

YES = 1 
) BOOL; 

typedef enum 

LO { 

SUCCEED = 0, 

FAIL = 1 

> SUCCESS ; . 

.5 =ifndef abs 

^define abs(x) (((x) < 0) ? -(x) : (x) ) 

^endif 

#ifndef max 

#define max(x,y) ( ( (x) < (y) ) ? (y) : (x)) 
0 =endif 

#ifndef min 

#define min{x,y) ( ( (x) < (y) ) ? (x) : (y) ) 
5=endif 

5 =define FOREVER while (YES) 



??def ine EOS ( ' \0 ' ) 
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WE CLAIM: 

1. A method for electronically compressing 
handwritten data expressible as a sequence of data points, the 
method comprising the steps of: 

capturing a parametric representation of data 
handwritten on a digitizer tablet; 

determining local extrema and saving representations 
of said local extrema as first data points; 

determining stroke endpoints and saving ' 
representations of said stroke endpoints as second data points; 
and 

calculating curvature locally about selected data 
points between said second data points and saving only selected 
representations of said selected data points based on 
curvature. 



2. The method according to claim 1 wherein said 

parametric representation of the digitizer data is a sequence 

20 Of X and y coordinate values in a Cartesian coordinate system 
as a function of time. 

3. The method according to claim i wherein said 
local extrema determining step comprises selecting minimum x^ 

2o and y^ and maximum x^ and y^ according to the following 
relationships: 

minimum = x^ < x^^^j^ and x^ < x^.^^, 
maximum = x^ > x^^^^ and x^ > x^.^^, 
minimum = y^ < y.^^ and y^ < y^_^, 
maximum = y^ > y.^^ and y^^ > y._^, 
where 

Xi and y . over the range i = 0 to i = nsamps-2 is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

4. The method according to claim 1 wherein said 
stroke endpoint determining step comprises selecting a first 



30 
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'^endpoint of a stroke, Xq and y^. and selecting a last endpoint 

of a stroke, x„^__„^ , and v 

nsamps-1 J^nsamps-1' 

where 

and over the range i = 0 to i = nsamps-l is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

5. The method according to claim 1 wherein said 
curvature determining step comprises calculating an angle 
between adjacent data points according to the relationship: 



~ ^i+d " ^i-d' 



where 



= - 360 degrees if > 180 degrees, 
= + 360 degrees if < -180 degrees, 
^i+d = a^ctan (Xi^^ - y^^^ - y . ) , 
^i-d = arctan (x^ - x^.^, y. - y . , 
d is an integer no larger than 4, 

x^ and y^ over the range i = d to i = nsamps-l-d is 
the discrete representation of a stroke, and 
nsamps is the number of points. 

6. The method according to claim 5 wherein said 
curvature angle calculating step further ccr.prises discarding 
data points with a curvature angle less than 5 degrees or 
greater than 25 degrees. 

7. The method according to claim 5 wherein said 
curvature angle calculating step further comprises discarding 
data points with a curvature angle less than 25 degrees or 
greater than 50 degrees. 

8. The method according to claim 1 further 
including a preprocessing step of discarding duplicate data 
points according to the following relationships: 

if >«i = x^.T and y^^ = y^_^, 
then discard x^^ and y^ ; 
if ^i - ^i-i and y. - y._^, 
then save x^ and y^ ; 
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'Where 

Xi and over the range i = 1 to i = nsamps-l is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

9. The method according to claim 1 further 
including a preprocessing step of discarding data points 
according to the following procedure: 

set j = 1, 

calculate d = { (Xq - Xj)2 + (y^ _ y.)2^h^ 

if d < 5, then discard point corresponding to index 

j, set j = j + 1, and recalculate d; 

If d > 5, then shift remaining points down so that 

the point (Xj, yj) corresponds to (x^, y^) , and quit the 

calculation, 

where 

x^ and y . over the range i = o' to i = nsamps-l is the ' 
discrete representation of a stroke, and 
nsamps is the number of points. 

10. The method according to claim 1 further 
including a preprocessing step of discarding data points 
according to the following procedure: 

set j = 2, 

calculate d = ( (x -y \2. 

* ^ -^nsamps-l ^nsamps-j) 

^ynsamps-i'^nsamps- j ) ^ 5 " ' 
If d < 5, then discard point corresponding to index 

j, set j = j + 1, and recalculate d; 

If d > 5, then shift remaining points down so that 

the point (Xnsan,ps-j» ynsamos-j) corresponds to (x_^^„„^ 

yn3air,ps-2^ quit the calculation. 

11. The method according to claim 1 further 
including a preprocessing step of dividing data points between 
said first data points into groups of three adjacent points and 
calculating the median value of said group according to the 
following relationships: 

X. = median value of (Xj_^, Xj , Xj^^), and 
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'• yj = median value of (Yj.^, yj , Yj+i). 

where 

and y^ over the range i = 1 to i = nsamps-2 is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

12. The method according to claim i further 
including a preprocessing step of dividing data points between 
said first data points into groups of three adjacent data 
points and calculating a representative value for said group 
according to the following relationships: 

Xj = 0.25Xj_;^ + O.SXj + 0.25Xj + ;L'^"d 
yj = 0.25yj_^ + o.5yj + 0.25yj^^. 

where 

x^ and y^ over the range i = 1 to i = nsamps-2 is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

13. An apparatus for electronically compressing 
handwritten data expressible as a sequence of data points, the 
apparatus comprising: 

means for capturing a parametric representation of 
data handwritten on a digitizer tablet; 

means for determining local extrema and saving 
representations of said local extrema as first data points; 

means for determining stroke endpoints and saving 
representations of said stroke endpoints as second data points; 
and 

means for calculating curvature locally about 
selected data points between said second data points and saving 
only selected representations of said selected data points 
based on curvature. 

14. The apparatus according to claim 13 wherein said 
parametric repress ation of the digitizer data is a sequence 
of x and y coordir. - values in a Cartesian coordinate system 
as a function of t:..-e. 
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apparatus according to claim 13 wherein said 
local extreina determining means comprises means for selecting 
minimum and and means for selecting maximum x^ and y. 
according to the following relationships: ^ 

minimum = x^ < x^^^ and x^ < x^.^^, 

maximum = x^ > x^^^ and x^ > x^.^, 

minimum = y^ < y.^^ and y^ < y^.^, 

maximum = y^ > y.^^ and y^ > y._ 

where 

Xi and yi over the range i = o to i » nsamps-2 is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

16. The apparatus according to claim 13 wherein said 
stroke endpoint determining means comprises means for selecting 
a first endpoint of a stroke, Xq and y^,, and means for 
selecting a last endpoint of a stroke, x.^^^^^ . and y 

where nsamps-l ^nsamps-l' 

Xi and y^ over the range i = 0 to i = nsamps-l is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

17. The apparatus according to claim 13 wherein said 
curvature determining means comprises means for calculating an 
angle between adjacent data points according to the 
relationship: 



where 



~ ^i+d ~ ^i-d' 



Ci - - 360 degrees if > iso degrees, 
Ci = + 360 degrees if < -180 degrees, 
ai^d = arctan (x^^^ - x^, y.^^ - y . ) , 
ai.d = arctan (x^ - x^.^. Yi " Yi.^) . 
d IS an integer no larger than 4, 

Xi and y . over the range i = d to i = nsamps-l-d is 
the discrete representation of a stroke, and 
nsamps is the number of points. 
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18. The apparatus according to claim 17 wherein said 
curvature angle calculating means further comprises means for 
discarding data points with a curvature angle less than 5 
degrees or greater than 25 degrees. 

19. The apparatus according to claim 17 wherein said 
curvature angle calculating means further comprises means for 
discarding data points with a curvature angle less than 25 
degrees or greater than 50 degrees. 



20. The apparatus according to claim 13 further 
including a preprocessing means for discarding duplicate data 
points according to the following relationships: 

if ^i = ^i-i and Vi = Yi.^, 
then discard and y^; 
if • x^.T and y^ * y^_^, 
then save x^ and y^; 

where 

Xi and y^ over the range i = i to i = nsamps-l is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

21. The apparatus according to claim 13 further 
including a preprocessing means for discarding data points 
according to the following procedure: 

set j = 1, - 

calculate d = { (Xq - Xj)2 + (y^ _ y.^l^h^ 

if d < 5, then discard point corresponding to index 
j, set j = j + 1^ and recalculate d; 

if d > 5, then shift remaining points down so that 

rhe point (Xj, y.) corresponds to (x,, y,) , and quit the 

calculation, 
where 

and y over the range i = 0 to i ^- nsamps-l is the 
discrete represenr^rion of a stroke, and 

nsamps is the number of points. 
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" 22. The apparatus according to claim 13 further 
including a preprocessing means for discarding data points 
according to the following procedure; 
set j =2, 

calculate d = < (x . ^ 2 . 

* ^*nsanips-l ^nsamps-j' 

^ynsamps-l'^nsamps-j ) 
if d < 5, then discard point corresponding to index 

j, set j = j + 1, and recalculate d; 

if d > 5, then shift remaining points down so that 

the point (x„s3„ps_j, ynsamps-j) corresponds to (Xj,samps-2' 

i'nsamps-2) '^^^^ calculation. 

23. The apparatus according to claim 13 further 
including a preprocessing means for dividing data points 
between said first data points into groups of three adjacent 
points and means for calculating the median value of said group 
according to the following relationships: 

Xj = median value of {Xj_;^, Xj , Xj^^^), and 
yj = median value of {y^.^, y . , y.^^} . 

where 

x^ and y^ over the range i = 1 to i = nsamps-2 is the 
discrete representation of a stroke, and 
nsamps is the number of points. 

24. The apparatus according to claim 13 further 
including a preprocessing means for dividing data points 
between said first data points into groups of three adjacent 
data points and means for calculating . a representative value 
for said group according to the following relationships: 

Xj = 0.25Xj_^ + o.5Xj + 0.25Xj^^, and 
yj = 0.25yj_^ + 0.5yj + 0.25yj^^. 

where 

and y . over the range i = 1 to i = nsamps-2 is the 
discrete representation of a stroke, and 
nsamps is the number of points. 



wo 94/03853 



PCT/US93/06883 



1/15 




NORMALIZATION/ 
PREPROCESSING 



REMOVE DUPLICATE/ 
CLOSE POINTS 



STAGE ONE 
FEATURE BASED 
COMPRESSION 



STAGE TWO 
MATHEMATICAL 
COMPRESSION 



-130 



-140 



-150 



■160 



vec ^ DOES 
Y-iS^^STAGE 2 IMPROVE 
COMPRESSION 



■170 



190- 



SAVE DATA WITH 
STAGE 1 AND STAGE 2 
COMPRESSION 



SAVE DATA WITH 
STAGE 1 
COMPRESSION ONLY 



-180 



F/G. 7 



wo 94/03853 



PCr/US93/06883 



2/15 



RANGE i = 0 TO i = NSAMPS - 2 



•210 



is1 



■220 



MINIMUM = Xi < Xm and Xj < Xm 
MAXIMUM = Xj > Xi+1 AND Xj > Xm 

MINIMUM = Yj < Ym and Y| < Ym 
MAXIMUM = Yj > Yi+1 AND Yj > Ym 



■230 




NO^ 

-^^^^ i = NSAMPS - 1 

9 



END 



FIG. 2 



wo 94/03853 



3/15 



PCT/l)S93/06883 




FIG. 3 




FIG. 4 



5/15 



PCT/US93/06883 



X- Y 

RANGE i = 0 TO i = NSAMPS - 1 



I 



•510 



i» 1 
XsAV - Xo 



-520 




DISCARD 
Xi. Yj DATA POINT 



•550 




END 



>-540 


1 T 
XSAV = Xj 




YSAV = Yj 



■545 



FIG. 5 



wo 94/03853 PCT/US93/06883 

6/15 



640 




635 



FIG. 6 



wo 94/03853 



PCr/US93/06883 



7/15 



X- Y- 

RANGE i = 0 TO i = NSAMPS - 1 



710 



i«0 



J = 1 



•715 



•720 



CALCULATE 
d = [(Xi-Xj)2 + (Yi-Yj)2]^/2 




DISCARD DATA POINT 
CORRESPONDING 
TOXj,Yj 



J = J + 1 



■760 



•730 



750 



SHIFT POINTS 
SOTHATXj.Yj 
CORRESPONDS TO 
Xi.Yi 



•770 



END 



FIG. 7 



wo 94/03853 PCr/US93/06883 

8/15 



RANGE i = 0 TO i = NSAMPS - 2 
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CALCULATE 
a = MEDIAN VALUE OF 

{Xj_i , Xj, Xj+i } 
b = MEDIAN VALUE OF 
{Yi-i.Yj,YM} 
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RANGE i = 0 TO i = NSAMPS - 2 



i» 1 
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Xi = 0.25 Xi_i + 0.5 Xi + 0.25 Xj+i 
Yi = 0.25 Yi_i+ 0.5 Yi + 0.25 Yui 
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FIG. 10a 
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FIG. 11a 
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FIG. 12a 
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FIG. 13a 
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FIG. 13b 



FIG. 13c 
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FIG. 14e 
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